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A layered superconductor with a full pairing energy gap can be driven into a nodal superconducting 
(SC) state by inter-layer pairing when the SC state becomes more quasi-3D. We propose that this 
mechanism is responsible for the observed nodal behavior in a class of iron-based SCs. We show 
that the intra- and inter-layer pairings generally compete and the gap nodes develop on one of the 
hole Fermi surface pockets as they become larger in the iron-pnictides. Our results provide a natural 
explanation of the c-axis gap modulations and gap nodes observed by angle resolved photoemission 
spectroscopy. Moreover, we predict that an anti-correlated c-axis gap modulations on the hole and 
electron pockets should be observable in the 5' ± -wave pairing state. 



For the iron-based superconductors [2 [2] with a com- 
plicated band structure, the symmetry of the order pa- 
rameter in the superconducting (SC) state [5] remains 
elusive. The 5 ± -wave pairing symmetry, predicted by 
both strong [3HZ] and weak coupling theories [5HTU] based 
on the magnetic origin, is a promising candidate and has 
been supported by many experimental results [6l [TTHT3] . 
However, it has also been seriously challenged by the exis- 
tence of gapless excitations or nodal behavior observed in 
some iron-based superconductors [HHI1|, in particular, 
BaFe 2 As2- x P x [20-23, 25 where some of the As atoms 
are replaced by the P atoms. A possible explanation of 
the nodal behavior has been suggested by the weak cou- 
pling approaches, such as the functional renormalization 
group (FRG) technique [27J HH] and random phase ap- 
proximations (RPA) [25] . These calculations suggest that 
gap nodes can develop on the electron pockets when the 
detailed nesting properties vary among the hole pockets 
located at the T and M points, and the electron pockets 
located at the X point of the unfolded Brillouin zone. 
When the size of the hole pocket at M point decreases, 
the SC gap on the electron pockets becomes increasingly 
anisotropic and eventually gap nodes emerge. The re- 
duction of the M hole pocket can be achieved by either 
increasing electron doping or by tuning the pnictogen 
height through replacing As by P [23"I l2l?] . 

Recently, nodes in the gap function dispersion along 
in the c-axis (c-axis nodal lines) have been observed di- 
rectly by ARPES in BaFe 2 As 1 . 7 P .3 [30 J. However, the 
weak coupling theories cannot explain the observed nodal 
behavior for the following reasons. First, the observed 
nodes are on the hole pockets, not on the electron pock- 
ets. Second, in contrast to LDA calculations [25], the 
P substitution in these materials does not push the hole- 
like band near M to sink below the Fermi surface [2"Tll25| . 
Instead, as the substitution increases, the M hole pocket 
and the X electron pockets barely change while the hole 
pockets near the T point at the zone center (k z = 0), 



which have large c-axis dispersions, accommodate the 
additional holes. As a result, with increasing P sub- 
stitution, one of the two T hole pockets grows increas- 
ingly larger. The size of this hole pocket at the Z point 
(k z — 7r) can even be larger than the size of the largest 
hole pocket in ^^62^2, the most hole-doped iron-based 
superconductors known today. These properties point 
to a non-rigid band picture under the "iso-valent" dop- 
ing and completely violate the assumption of the band 
structure taken in the above weak coupling theories. 

In this Letter, we suggest that the observed nodal be- 
havior originates from the inter-layer pairing and the re- 
duction of the intra-layer SC pairing gap due to the in- 
crease of the size of the hole pockets. This proposal con- 
sistently explains the c-axis modulation of the SC gaps 
observed in optimally hole-doped Bai- x K x Fe2As2[3l, 
132] and the nodal behaviors in BaFe2As2- x P x ]JE}\- It 
suggests that the gap modulations along the c-axis are 
directly related to the c-axis band dispersion. Our cal- 
culation also reveals that the inter-layer SC pairing, in 
general, competes with the intra-layer SC pairing in 
these quasi- two dimensional materials, a possible rea- 
son why the highest T c is not achieved in the 122-family 
(AFe 2 As 2 ) but in the 1111-family (AOFeAs) of the iron 
pnictides |341 135] since the former is more three dimen- 
sional [2 [331 EZ] ■ Moreover, we predict that the S' ± -wave 
pairing symmetry should result in an anti-correlation of 
the SC gap values between the hole and electron pockets 
along the c-axis as a function of c-axis momentum. This 
property, if observed, can serve as a direct experimental 
evidence for the S -wave pairing symmetry. 

Model We construct a three-orbital model which in- 
cludes the d xz , d yz and d z ? orbitals to study the physics. 
Experimentally, the large c-axis dispersion is only ob- 
served in one of the hole pockets near the T point which 
is mainly composed of d xz ^ yz orbitals [36] • The increase 
of the c-axis dispersion upon P doping is mainly due to 
the increase of the mixture of the d z i orbital into this 



2 



hole pocket [35]. This has been shown by both polar- 
ized ARPES experiments [25] and numerical calculations 
[291 EH]. The ARPES experiments [25j [30] show that 
the hole pocket with the large c-axis dispersion has even 
symmetry with respect to the reflection of the T — M 
mirror plane and, with increasing P doping, the band 
mainly attributed to the d z 2 orbital moves closer and 
closer to the Fermi energy so that the weight of d z 2 on 
the Fermi surface increases. This picture is consistent 
with the symmetry analysis since the d z 2 orbital is also 
symmetric with respect to the reflection of the T — M 
mirror plane. The model we construct captures all the 
above essential experimental results and can still achieve 
high analytical tractability. We will also show that the 
results for the c-axis properties derived from this model 
are rather generic. 




FIG. 1: 3-D Fermi surfaces and Contour plot of Fermi sur- 
faces in the extended three-band model Eq.|TJ (fc z = for the 
read line and k z = n for the green line). 

The Hamiltonian of our model includes two parts H = 
H t + H], where H t is the kinetic energy and Hi is the 
pairing interaction. H t is given by 

H t = Y2 £ kal3dl aa dui3a, (1) 
k.a/3<T 

where a, (3 = 1,2,3 label the 3d electrons in the d xz , 
d yz and d z 2 orbital respectively. The form of the kinetic 
energy and the hopping parameters are constructed by 
extending the two-orbital model [42] to achieve the c- 
axis dispersion that matches well the experimental re- 
sults. Their explicit forms are 

£k,n = — 2ti cos k x — 2i2 cos k y — 4^3 cos k x cos k y — /i, 
£k,22 = —2^2 cos k x — 2ii cos k y — 4^3 cos k x cos ky — /i, 
£k,33 = 2i z (cos k x + cos k y ) - fi 3 , 
£k,i2 = —4*4 sin fc^ sin fey, 

£k,i3 = —txz (1 - cosfc z ) (cosfer + cos k y ) sinfe^,, 
£k,23 = txz (1 - cosfc z ) (cos + cos fey) sinfe y . 

The c-axis dispersion is induced by the coupling between 
the d xz , yz and d z 2 orbitals. The later is taken to be be- 
low the Fermi level. In the explicit form of the coupling 



between these orbitals, £k,i3, we have taken into account 
both the lattice symmetry requirements and the exper- 
imental observations. By taking the following parame- 
ters t\ = — 1,< 2 = 1-55, t 3 = i 4 = — 0.85, t z — l,t xz — 
0.8, n — 1.77, /13 = 5, we find that the model describes 
well the 3-dimensional Fermi surfaces measured experi- 
mentally (25[ EH] as shown in Fig. [I] 

Hj, the SC pairing interaction, generally includes the 
following terms, 

HintraJ = ~Vl ^ (f>^ Vk^4at^-kc4 C '- k ' c 4^ k ' a T 
kk'.Q=l,2 

Hi n ter,2 = ~V% ^ 0fc 0k' 4at^-kc4^-k'o4-^k'a!t 
kk',a=l,2 

#mter,3 = ~V 3 ^ ^kVk^Sf^-k.S^-^S^k'.St 
kk' 

Hi n ter,A = — Vj ^ ^Vk' (^k.lt^-k.ll^-k',34.^k',3t 
kk' 

+rf k; 2t rf -k,24.^-k',34.rfk',3t + h.C.). 

The first term Hi ntr a,i describes the intra-layer pairing 
while the rest three terms describe the different inter- 
layer pairing interactions. H inter 2 accounts for the inter- 
layer pairing interaction between the d xz and d yz or- 
bitals, as does Hi nter ,3 for the d z 2 orbital. The last 
term describes the inter-layer interaction between the 
d z 2 and d xzyz pairs. We take <j>^ to be the intra- 
layer S ,± -wave paring function (f>^ = cosfc^cosfcj, and 

= (cosk x -t- cos ky ) cos k z , and k 3,4 ' ) = cosfe 2 . These 
choices rely on the assumption that the SC pairing in 
iron-based superconductors is rather short-ranged. The 
form of (j)^ has been proposed in the models based on 
local magnetic exchange couplings [IHI] and it has been 
shown that the form factor is consistent with current ex- 
perimental results [TT]. The form of 4>^ has been pro- 
posed in [35], which can be obtained from the existence 
of AFM exchange couplings between the layers E2 HO] • 
The form of </>( 3 ' 4 ) and the corresponding V34 pairing in- 
teractions can be understood as the inter-layer pairing 
is between two adjacent layers and the pairing symme- 
try is s-wave. V3 describes the inter-layer pairing for the 
d z 2 orbital. The V4 term, which describes the coupling 
between two inter-layer pairings of two different orbitals, 
can be understood in the following way. Since the c-axis 
hopping term in H t describes the hopping between two 
adjacent layers and the d z 2 is below the Fermi level, the 

second order perturbation through such hopping would 

t 2 

generically produce V4 oc 

In the self-consistent mean-field theory for the SC 
state, Hj becomes 

Hbcs = £ A a (k)4 a /_ w + h.c, (2) 

ka 

where A x (k) = + a£ 2) ^ 2) + A< 4) k 4) , A 2 (k) = 
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A<M 2) + ArC ; A 3 (k) = A^>^ 



(A[ 4) + A^ 4) )0k } - Here aL™ } is denned by 



A<"> = V n 



k 



(3) 



Results Before we present a full numerical solution for 
the above Hamiltonian, we first discuss the simple phys- 
ical picture for the generation of the nodal points in the 
gap function on the hole pocket. In the above model, if 
we consider the general intra-orbital pairing form of the 
d X z,yz orbitals, the /c z -dependent SC gap can be written 
as 

A(k) = Ao [cos k x cos k y + 5 Z (A + cos k x + cos k y ) cos k z ] , 

(4) 

where the first term represents the S pairing and the 
second term represents the inter-layer pairing with the 
S'-wave pairing symmetry. In the first order approxima- 
tion, this intra-orbital pairing roughly determines the SC 
gap since it dominates as we will show later. This form 
indicates that the inter-layer pairing is between the two 
neighboring layers. The gap zero points develop as S z 
increases. As shown in Fig. [2j when S z reaches a certain 
value, the contour of the gap zeroes will cross the Fermi 
surface at the points near k z = ±7r, which leads to the 
nodal behavior. 




FIG. 2: A cartoon plot for the gap nodal picture in the 
momentum space (k x ,0,k z ). The blue lines are the contour 
lines for zero gap value of Eq. Q with 8 Z = 0.4 and A = 0. 
FSi and FS2 are defined in Fig. ITT The intercept between the 
blue and the red lines produce the nodes in the hole pocket. 



Second, we discuss two important, general results ob- 
tained from our model, which are independent of the de- 
tailed pairing interaction parameters Vi in Hj. One of 
these is that the inter-layer pairing always competes with 
the intra-layer pairing. To demonstrate this more clearly, 
we switch off V3 and V4 and perform a self-consistent solu- 
tion with V\ and V 2 . The SC pairing gaps as a function 
of the interacting parameters are shown in Fig|3] It is 
very clear that the intra-layer pairing reduces while the 
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FIG. 3: Gaps versus pairing interactions. Interaction pa- 
rameters in (a) are defined as V2 = 3, V3 = V4 = and (b) 
Vi = 3, V 3 = Va = 0. Ai n) is defined in Eq. jil. 




FIG. 4: Gaps versus pairing interaction V2 with Vi = V3 

V4, = 3. 



inter-layer pairing increases and vice versa. If we turn 
on all of the inter-layer pairing interactions, the results 
are rather similar: while the different inter-layer pair- 
ings can increase simultaneously, the intra-layer pairing 
gaps always decrease as the inter-layer ones increase. A 
typical result is presented in Fig(4] This result quali- 
tatively suggests that a more quasi two-dimensional SC 
state is likely better for achieving a higher T c since the 
intra-layer pairing would dominate. So far the highest 
T c in the iron-based superconductors is achieved in the 
1111-family. The highest T c in the 122 family is about 8 
degrees lower than the one in the 1111-family [2]. Com- 
paring to the 122 family, the 1111 family is much more 
two-dimensional with much less dispersion along the c- 
axis. 

In a multi-orbital model, the relation between the SC 
pairing parameters and the energy gap in the low energy 
single particle excitations can be complicated. In order to 
show that the SC state truly develops nodes, we have to 
calculate the energy dispersion of the Bogoliubov quasi- 
particles at the Fermi surfaces. In FigjEJ we plot the dis- 
persion of Bogoliubov quasiparticles along the c-axis with 
the parameters given by V% = 3, Vi = 3.35, V3 = V4, = 3. 
There are two important results. One is that the true 
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nodes can easily develop on the hole pocket. The other 
is that the gap values of the quasiparticles on the hole 
pocket and electron pocket are anti-correlated along the 
c-axis: the gap value on the hole pocket is larger at k z = 
than at k z = tt while on the electron pocket, it is smaller 
at k z — than at k z = tt. This anti-correlation is a com- 
bined result of the inter-layer pairing and the S^-wave 
symmetry for the intra-layer SC pairing order parame- 
ter which is proportional to cosk x cosk y and changes sign 
between the hole and the electron pocket. This result 
holds for most of the parameter regions we have inves- 
tigated. Therefore, this inter-layer pairing induced anti- 
correlation suggests that ARPES can provide a direct 
test of the possible 5' ± -wave pairing gap functions in the 
iron-pnicitide superconductors. Of course, to detect it, 
a high energy-resolution in the ARPES experiments has 
to be achieved since the c-axis dispersion and the gap 
modulation on the electron pockets are not large. 




FIG. 5: The energy dispersion of the Bogoliubov quasi- 
particles at the Fermi surfaces shown in Fig. [I] Here 



Vi = 3, V 2 = 3.35, V 3 = 14 = 3 with A[ 1} 2 = 0.078, A\% = 



,(2) 



0.057, Af ] = 0.007, A$ = 0.016, = 0.007. 

Finally, we note that several previous thermal con- 
ductivity measurements have suggested that the nodes 
should be on the electron pockets [TS]. A key argument 
given in [T5] is that the quaisparticles on the hole pockets 
have much lower velocity and heavier mass than those on 
the electron pockets so that the in-plane Fermi velocity 
Vf of the hole pockets is too small to explain the observed 
residual thermal conductivity. However, this statement 
is only partially true. There are three hole pockets cen- 
tered around the folded Brillouin zone center. The Fermi 
velocity on one of the hole pockets is in fact comparable 
to that on the electron pockets. ARPES results [33] show 
that the former is even slightly larger than the later. This 
hole pocket, whose orbital character is even with respect 
to the r — M mirror plane, is exactly the pocket that 
carries the large c-axis dispersion. Therefore, the pre- 
vious thermal conductivity measurements are consistent 
with our results for the existence of gap nodes on the hole 
pocket. 



,(4) 



(4) 



Conclusion We have constructed a model to show how 
nodes in the single-particle excitations can emerge in 
iron-based superconductors. The development of the 
nodes are due to the combined effects of the increase 
in the hole pocket size which reduces the SC gap from 
intra-layer pairing and the presence of the inter-layer SC 
pairing. This study consistently explains the experimen- 
tal observations of the c-axis gap variation in optimally 
hole-doped Bai^ x K x Fe 2 As2 [22 E2] and the nodal be- 
haviors in BaFe 2 As2- x Px ISO]- We also demonstrated 
that the inter-layer and intra-layer pairing generally com- 
petes with each other, and suggested a direct experi- 
mental test of the 5' ± -wave pairing symmetry through 
the anti-correlation of the gap modulations on the hole 
and electron pockets that can be measured by ARPES. 
We believe that our results can also explain the observed 
nodal behaviors in other materials such as -£^62^432 and 
AFe2- x Ru x As2- A concrete test of our model will be 
whether the gap functions observed in these materials 
obey Eq. ^ as the leading contribution to the quasi-3D 
SC pairing gap function. 
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